IHES/P/95/99 



Testing for gravitationally preferred directions 
using the lunar orbit 



Thibault Damour 

Institut des Hautes Etudes Scientifiques, 91440 Bures sur Yvette, France 
and Departement d'Astrophysique Relativiste et de Cosmologie, Observatoire de Paris, 
Centre National de la Recherche Scientifique, 92195 Meudon CEDEX, France 

David VokrouhlickyQ 

Observatoire de la Cote d'Azur, Departement CERGA, Avenue N. Copernic, 

06130 Grasse, France 
(February 7, 2008) 



Abstract 

As gravity is a long-range force, it is a priori conceivable that the Uni- 
verse's global matter distribution select a preferred rest frame for local grav- 
itational physics. At the post-Newtonian approximation, the phenomenology 
of preferred-frame effects is described by two parameters, a.% and «2, the sec- 
ond of which is already very tightly constrained. Confirming previous sugges- 
tions, we show through a detailed Hill-Brown type calculation of a perturbed 
lunar orbit that lunar laser ranging data have the potential of constraining 
a,\ at the 10 level. It is found that certain retrograde planar orbits exhibit 
a resonant sensitivity to external perturbations linked to a fixed direction in 
space. The lunar orbit being quite far from such a resonance exhibits no sig- 
nificant enhancement due to solar tides. Our Hill-Brown analysis is extended 
to the perturbation linked to a possible differential acceleration toward the 
galactic center. It is, however, argued that there are strong a priori theoretical 
constraints on the conceivable magnitude of such an effect. 



1 



Typeset using REVTeX 



I. INTRODUCTION 



It has been recognized since many years that the lunar motion provides a superb testing 
ground for relativistic gravity In particular, the Lunar Laser Ranging (LLR) experi- 

ment has allowed one to get a very high precision test of the equivalence principle, as well as 
a 1 % test of the einsteinian spin-orbit coupling |7|||. However, it has been recently pointed 
out that the lowest-order perturbation analyses that have been commonly used [£|,[|,|| to 
derive theoretical estimates of (null or non-null) relativistic effects are insufficiently accurate 
in view of the importance of solar tidal effects [[Hj. Motivated by the results of Ref. UTOj , 
we presented in Ref. [JTTJ] a full-fledged Hill-Brown theory of the lunar orbit perturbation 
due to an hypothetical violation of the equivalence principle. We found that the interaction 
with the quadrupolar tide amplified the results of lowest-order perturbation analyses by a 
very significant factor: 60 % increase of the naive first-order calculation, or 40 % increase 
of the improved first-order calculations allowing for perigee motion. Such results raise the 
question of whether similar amplification factors affect other (null or non-null) relativistic 
effects in the lunar motion. To address this question it is convenient to use the parametrized 
post-Newtonian (PPN) framework (see e.g. |9[]) in which possible deviations from general 
relativity in the weak- field regime are described by some parameters, (3 — 1, 7 — 1, a%, oli, 
etc., which vanish in Einstein's theory. 

In the case of the effects linked to the Eddington post-Newtonian parameters (3 and 7 
Ref. |5| has indeed shown that tidal effects are numerically important. An observationally 
oriented discussion of the influence of the tidal deformation on the main effects linked to (3 
and 7 is contained in Ref. fl0|| , while, as we said above, the tidal amplification of equivalence- 



principle- violation effects was discussed in Refs. [10] and, in more detail, in [] 1 



In the present paper, we study the influence of the tidal deformation of the lunar orbit 
on the preferred- frame effects linked to the parametrized post-Newtonian parameter a.\. 
We shall also discuss the effect of an hypothetical violation of the equivalence principle of 
galactic origin (noting, however, that there are strong a priori theoretical constraints on the 
magnitude of such a violation). 



II. PREFERRED FRAME EFFECTS 

As gravity is a long-range force, one might a priori expect the Universe's global matter 
distribution to select a preferred rest frame for local gravitational physics. As shown in 
T2] , p!3| j9[] , all preferred frame effects in the first post-Newtonian limit are phenomenologically 



describable by only two parameters, «i and 0.2- These parameters are associated with the 
following terms in the Lagrangian describing the gravitational dynamics of A-body systems 

(a,b 1 V) 

a x x . Gm A m B / o\ / Q -, x 

= - („ AB .v») (n AB .v%)] . (2.1b) 

4 A+B r ABC 
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Here, represents the velocity of a given body with respect to the gravitationally preferred 
frame and n^s = (r^ — r#) / r ab- Many (though not all) of the observable effects linked to 
a,\ and a.2 depend on the choice of the gravitationally preferred frame. We shall follow the 
standard assumption [[J that the latter frame, being of cosmological origin, can be (at least 
approximately) identified with the rest frame of the cosmic microwave background. This 
means that the center of mass of the solar system has the velocity w with respect to the 
preferred frame of rest, with |w| ~ 370 ± 10 km/s in the direction (a, 5) = (168°, —7°) | ]14| . 

It has been shown in Ref. [15| that the close alignment of the Sun's spin axis with the 
solar system angular momentum yields an extremely tight bound on a^- \o.^\ < 3.9 x 10~' 
(90 % C.L.). This limit on is much stronger than the existing limits on the other post- 
Newtonian parameters (3, 7 and ot\. We shall therefore neglect all «2 effects in this work. 
Concerning a±, combined orbital data on the planetary system yield |IB[ 



e*i = (2.1 ±3.1) x 10" 4 (90%C.L.), (2.2) 
while binary pulsar data yield comparable or better limits ||17||. More precisely, PSR 1855 



+ 09 data yield < 5.0 x 10~ 4 (90%C.L.) 0], while a recent analysis of PSR J2317 + 
1439 data |T§ yield 

|ai| < 1.7 x 10" 4 (90%C.L.) . (2.3) 

The fact that the observational limits on the ot\ parameter are only a factor ten better 
than the present limits on the (more conservative) Eddington post-Newtonian parameters 



(3 — 1 and 7 — 1 stimulated recently Damour and Esposito-Farese JT9| to propose several 
experiments for improving them. Concerning their proposal to use artificial satellite motions, 
it has been recognized that the currently best laser tracked satellite LAGEOS cannot yield 
a better constraint on a± because of badly modeled non-gravitational forces [M. Another 
possibility mentioned in Ref. jlJj (and first pointed out in ||) concerns the lunar motion and 
suggests that LLR data might yield an interesting new limit of the ot\ parameter. However, 



Refs. H and |L9| used only first-order perturbation theory to estimate the ai-effects in 
the lunar motion. In view of these facts (and the experience of the strong coupling with 
the solar tides mentioned in Sec. I), we decided to reassess the quantitative value of lunar 
data for constraining ot\ by building an accurate Hill-Brown theory of the preferred-frame 
perturbations of the lunar orbit. 



We thus consider the 3-body Earth- Moon-Sun system (keeping the notation of Ref. |TT 
in particular we use the labels 1 = Moon, 2 = Earth, 3 = Sun). Generally all these three 
bodies contribute to the sum in Eq. ( |2.1a| ), however, we shall restrict ourselves to the "direct" 
preferred frame effects with the subscripts A and B spanning only 1 (Moon) and 2 (Earth). 
It is easy (though not trivial) to verify that the "tidal" preferred frame effects, involving 
the subscript 3 (Sun) in Eq. (|2.1a|) , are several orders of magnitude smaller than the direct 
effects. Because of their observational irrelevance, we also omit from our discussion several 



terms in Eq. ( p.laj) which are equivalent to a nearly constant redefinition of the locally 



measured gravitational constant. The dominant preferred frame effects are then contained 
in the following three contributions to the lagrangian (we factorized the Earth-Moon reduced 
mass /Z12 = mim2 /mo, mo = m\ + m2, from the lagrangian) 
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= (vo . v) t (2 . 4b) 

fl£) = _|Gp (wVo) (24c) 

Here, r = xi — X2 is the geocentric lunar position vector and v = dr/dt its velocity, vo is 
the velocity of the Earth-Moon center of mass motion around the Sun, and X21 = X2 — X\, 
with the mass ratios X\ = mi /mo and X2 = mz/mo = 1 — X\. In the following section, 
we treat successively the perturbations of the lunar orbit associated with the three terms 



III. HILL-BROWN TREATMENT OF PREFERRED FRAME EFFECTS 

A. The method in brief 

Since the Hill-Brown approach to lunar motion represents a classic tool of celestial me- 
chanics treated with care in the literature (e.g. [21-p4]|) we outline its concept only very 



briefly, focusing mainly on the particularities of the method involved in the present study. 
We also refer the interested reader to Ref. [T1J for more details and used notation. 

Following Hill we start by considering the planar Earth-Moon-Sun 3-body problem with 
the Earth-Moon center-of-mass on a circular orbit around the Sun (the so called "Main 
Problem"). The near circular lunar motion is investigated in an Earth-centered coordinate 
system (X, Y) rotating with the angular velocity n' corresponding to the solar motion around 
the Earth-Moon center-of-mass (the Sun thus rests on the axis X). Apart from the Earth 
direct gravitational action, the quadrupolar piece of the solar (tidal) gravitational potential 
is also taken into account. The reduced lagrangian of the lunar motion then reads 



Lmw 

with 



= i (X 2 + Y 2 ) + n' (XY - YX) + ifem , (3.1) 



Rw = + ~n' 2 X 2 (3.2) 

VX 2 + Y 2 2 V ; 

(the overdot means d/dt). Although simplified, the theory entails the most important part 
of the solar tidal deformation of the lunar motion. Instead of the usual keplerian ellipse Hill 
chooses for the intermediary lunar orbit a periodic solution (with particular symmetries) 
of ( |3.1| ), the so called "variational curve". In the following, we shall investigate the forced 
perturbations of the variational orbit due to the additional lagrangian terms I|f1- 



1 We recall that the variational orbit is not a general solution of the system ( |3.1[) . Apart from 



4 



A very convenient parametrization of the Hill problem consists of replacing (X, Y) by 
the complex conjugated quantities (w,w) defined by {i is the complex unit) 

X + iY = aC(l +w) , (3.3a) 

X -iY = aC l (l + w) , (3.3b) 

C = e lT , (3.3c) 

r = (n — n')t + r , (3.3d) 

where n denotes the mean lunar motion around the Earth and where the r variable represents 
the mean geocentric angular separation of the Moon and the Sun. Following Ref. ||, the 
fiducial lunar semi- major axis a is defined by 

Gm ° K{m) , (3.4) 



(n — n') 2 a 3 
where 

K(m) = 1 + 2m + -m 2 . (3.5) 

and m = n'/(n — n') is Hill's expansion parameter. (In the actual case of the Moon, 
mMoon = 0.0808489375. . ..) The Lagrange dynamical equations of the variational motion 
then read 

L(w,w) = W m \\(w,w) , (3.6) 

where we denoted 

3 

L(w, w) = D 2 w + 2(1 + m)Dw + -K{m){w + w) , (3.7) 

with D = d/(idr) = (d/d(, and the Hill source terms 

3 

W Hm (w, w) = --m 2 (- 2 (l + w) + K(m)Q{w, w) . (3.8) 
The nonlinear source function Q(w,w) and its development in terms of (w, w) can be found 



for instance in Ref. 11 



When considering an extra perturbation of the lunar motion, such as (|2.4| ) in the case of 
preferred frame effects, we have to include an additional source function on the right hand 
side of Eq. ( |3.6|) given by 



the "forced" perturbations related to a new physical cause it admits also "free" perturbations 
covering the classical notion of the lunar orbit eccentricity and its perigee drift due to the solar 
action. For simplicity, we omit in this study a natural coupling of the two types of perturbations 
neglecting thus the (small) lunar eccentricity (and inclination) corrections to the preferred frame 
perturbations. 
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S (w, w; Dw, Dw) = D 



9 G 9 G 
dDw dw 



(3.9) 



in terms of the "generating function" G(w,w;Dw,Dw) = 2{m/n'a) 2 R. 

Our method of solution of the system (|3.6| ) to (|3~^ ) consists of consecutive iterations, 
where at each stage one constructs a particular right-hand side source based on the results 
of the previous iterations. Details can be found in Ref. || or Appendix B of Ref. ||11|| . Let 
us only point out that, contrary to the simpler case of the synodic lunar perturbations due 
to an hypothetic violation of the equivalence principle studied in Ref. JTl| , the generic form 
of the right hand side source term now reads 



W*(a) = W. a C a + W a C 



(3.10) 



where we allow for: (i) complex functions W- a and W a , and (ii) any real (non integer) values 
of the powers a. Inversion of the linear problem L(w*(a), w+(a)) = W*(a) (a ^ 0) has a 



simple form w*(a) 



W. 



1 

A a (m) 
1 

AJm) 



+ w a ( a with 



a 



a 



2 (1 + m) a + -k 



W a - -kW. 
2 



2(1 



m) a H — k 
' 2 



W. 



and 



AJm) 



a' 



a 2 + 3k - 4 (1 + m) 



(3.11a) 
(3.11b) 

(3.12) 



[The solution corresponding to a = is identical with that given in Eq. (2.52a) of Ref. fpTTf] .] 
Some values of the power a in (|3.12| ) may lead to a significant amplification of the effect 
due to the smallness of the corresponding denominator A a (m). Of particular interest for 
our present work is the case where a = 1 + m which yields the small denominator 



A 1+m (m) = -m 2 (l + mf 



(3.13) 



In the next section we shall see that it appears in the sidereal excitation of the lunar orbit. 

Because of the background motivation of our work, related to the LLR experiment, we 
are essentially interested in the perturbation of the radial geocentric distance of the Moon 
given by 



er(l +w){l + w) 



(3.14) 



Performing a variation of this quantity, keeping only linear terms in the perturbation, we 
obtain 



a 



1 + w\ 



■7,\V2 



1+W 



5w 



(3.15) 



for the searched perturbation in radial coordinate. Remembering that w = 0(m 2 ), to lowest 
order in the m parameter, the radial oscillation can be expressed by the simple formula: 
5r/a ~ (w + w)/2. 
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In the rest of this section we investigate the forced perturbations of the lunar varia- 
tional orbit related to the three preferred- frame lagrangian terms (|2.4j ). Finally, we note 
that albeit the iteration scheme mentioned previously is straightforward it represents a huge 
algebraic manipulation exercise. We thus employed the powerful dedicated algebraic com- 
puter system MINIMS developed by M. Moons from the University of Namur (Belgium) |25| 



to perform this task. The lowest two orders of the results have been, however, checked by 
hand computations. 

B. Potential r£} 



Firstly, we focus on the source term ( 2.4a| ). The corresponding generating function G 
reads 

G(w, Dw) = -iei (?) { [Dw + (1 + m) (1 + w)\ C 1+m e ^ + [Dw - (1 + m) (1 + w)} C {1+m) e lc ^} 

(3.16) 



where 

h = ^xj^ (1) , (3.17) 

2 cr \a' J m 

and where one must express r in terms of w through Eq. ( j3.14|) . Here, a' is the radius of the 
(circular) solar orbit in the Earth-Moon center-of-mass frame and vq = n'a! its (circular) 
velocity, C(~ 0.98) is the cosine of the ecliptic latitude of the unit vector w°, and is a 
longitude angle of w° measured from the lunar (and solar) position at time to corresponding 
to an arbitrary new-moon phase. For instance, if we choose the last new-moon phase in this 
century, occurring at MJD51 521.2, we obtain <fi = 267.2°. Inserting this expression into 
Q3.9p we obtain the source function, to be added to the right hand side of the Hill equation 



( p.6|) , in the following form 



S (w,w; Dw, Dw) = (~^j ^Dw 



P -i<t> p i<i> 
£l+m_ ^- (l+m) c 

1 + w 1 + w 



+ (1 + m) 



1 + w 



(3.18) 



Working out the iterative solution mentioned above one realizes that this perturbation 
yields a wide spectrum of radial and longitudinal oscillations of the lunar orbit (compare 
also with the less accurate solution in Ref. ||). However, a detailed analysis shows that only 
two of them are sufficiently amplified to give an observably interesting signal: (i) terms with 
frequency equal to the mean sidereal lunar motion n (having a period of about 27. d 32), and 
(ii) terms with frequency equal to n — 2n' (having a period of about 32. d 13). Both periods 
are evaluated for the lunar orbit. Hereafter we discuss properties of both of them starting 
with the sidereal terms. 

The perturbation series giving the sidereal-frequency radial oscillations of the lunar orbit 
reads 
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T = H^* {m) sin [n {t ~ to) ~ 0] ' (3 ' 19) 

with 

„m / n , 67 395 2 103007 , 

M = 1 - Y m + ~Y m — 38^ m 

3327349 , ^. 

+ ^304~ m + ° {m ) • (3 - 20) 

Table | gives the coefficients of the series (m) up to the ninth order. The second column 
of the table indicates the numerical contribution of the corresponding term to the total value 



of the series for the lunar orbit, i.e. m = my[ oon = 0.0808489375. . . |2"4* |. Two important 
features are to be noticed: (i) a significant contribution of the higher order corrections to 
the total value of the series SjV(myiooa)-> an d (ii) a geometric-like character of this series 
clearly pronounced after a few terms. The second property suggests the presence of a pole 
near the valueP]m cr ~ —0.18407. Taking a Pade approximant of the series S^(m) confirms 
the presence of a unique root of the denominator located at the value m cr = —0.18407 with 
an error of about 10~ 5 . 

The physical origin of this pole can be easily understood by using the following argument. 
When solving the problem by traditional first-order perturbation techniques (see e.g. H,|B|), 
one finds that the sidereal orbit oscillation induced by an external force linked to a fixed 
direction in space contains, in the denominator, the secular rate of the perigee longitude w. 
This agrees with the intuitive idea that a spatially "frozen" orbit (not moving its pericenter 
in fixed space) is "resonantly sensitive" to constant forces. As a check of this idea, we 
have multiplied the S^(m) series by w(m) = n |m 2 + ^[m 3 + • • • as given by Andoyer 

up to order m 9 in Ref. [p2fl . The resulting series, say S^(m) = AS^(m)zj(m)/(3nm 2 ) = 
1 — m + ^-m 2 + ■ • •, is much more tame, showing that the main characteristics of Sj^(m) 
are entailed in the factor [cb(m)] -1 . The difference of our more precise solution (|3.19|) 
with the previous ones PJT9| (when the latter are improved by using the full value of the 
perigee advance) is essentially contained in the value of the residual series S { £(m). We 

find that, in the lunar case, (ttimooh) — 0.956. In contrast with the case of equivalence- 
principle-violation effects, we see therefore that preferred-frame effects exhibit no significant 
enhancement genuinely linked to the tidal deformation of the lunar orbit. This results holds 
for the other effects discussed below and is basically attributable to the fact that the actual 
lunar orbit is quite different from the "spatially frozen" resonant orbit (while it is rather 
near the orbit resonant for solar-directed equivalence-principle- violation effects). 

The principal quantitative information of the previous analysis is given by the amplitude 
of the sidereal oscillation of the lunar orbit ( |3.19D p|. Employing the definition ( |3.4| ) of the 



2 Such a value corresponds to a retrograde orbit with a sidereal period (for Earth satellites) T = 
2vr/|n| = 82.4 days. 

3 Beware, however, that for the final determination of the a\ constraint through LLR data anal- 
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auxiliary lunar semimajor axis a [and a Pade approximant value of S^(m) : S'W(mMoon) 
0.5465] we have 



an Jw\v Q 
\5 e r\ ~ — X 21 C- 



3 " c 



2 



21 !/ 3 

m 5 



S^{m) a' ~ 4780 x Ql C [cm] , (3.21) 



where x = (1 + m 3 /m ) _1 . The current published accuracy of the lunar ranging measure- 
ments performed by the CERGA team is 14 millimeters. Recent technical improvements are 
giving a timing precision of about 6 millimeters (C. Veillet, private communication). If the 
latter precision level can be turned into an accuracy level, the result (|3.21|) suggests that 
the LLR data should soon be able to constrain ai\ at the 1 x 10~ 4 level or better (given the 
phase information and the presence of several cui-effects at different frequencies). 

The value m = of the Hill parameter is apparently another singularity of our "ranging 
formula" (|3.19|) . However, because we neglected the Earth quadrupole and the other higher 
multipoles of the Earth gravity field, we cannot extend our solution to near-Earth satellite 
orbits (m ~ 0). This regime has been thoroughly discussed in Ref. |19|]. In a first approxima- 



tion we can, however, match smoothly the solution of Ref. |19| , accounting basically for the 
Earth quadrupole, and the solution presented in the present study, accounting in detail for 
the third-body (Sun) perturbations, by adding the Earth quadrupole contribution to zu{m) 
after having factorized it as denominator of the series S^(m). [Actually, as Andoyer's series 
is not accurate enough to localize precisely the zero at m — m cr , we found better to add the 
quadrupole contribution to the denominator of a Pade approximant of S^(m).] Figure |l| 
shows the synthesis of the two effects. The arrow points the singularity m = m cr , while the 
two points M and L stand for the Moon and the artificial (retrograde) satellite LAGEOS, 
respectively (without taking into account the LAGEOS inclination). We can see that none 
of the two bodies (an equatorial satellite at the LAGEOS altitude or the Moon) is the best 
candidate for testing preferred frame effects, but that (as mentioned in |]19|| ) a high orbit 
artificial body with a period of about 30 hours optimizes the sensitivity to the ati parameter 
(among prograde orbits). On the other hand, one should be aware of the fact that the 
motion of artificial bodies is typically influenced by many non-gravitational forces, some of 
which are difficult to be predicted and/or carefully modeled. For instance, this is the reason 
why the LAGEOS satellite is currently less suitable for constraining the a\ parameter than 
the Moon [pOfl , which is a nearly perfect "drag free Earth-satellite". Therefore the lunar 
data stand out as a potentially important source for the study of preferred frame effects. 

The intricate interaction of the variational curve perturbations with the underlying tidal 
deformation leads also to a slowly convergent series for the perturbations at the (n — 2n')- 
frequency. The final result for the radial oscillations formula reads 

& 4~ = -7-^ ( ! } M sin [(n - 2ri) (t - t ) + 0] . (3.22) 
a 4m 



ysis, the particular phase 4> is very important. Moreover, we shall see that the preferred frame 
perturbations of the lunar orbit act also with several other frequencies and it is their combined 
influence which determines the full effect to be searched for in the LLR data. 
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with 

„, m , , 43 28867 , 468391 , 4N 

Sf) m = 1 m + m 2 m 3 + C m 4 . 3.23 

ai v ; 6 720 2160 v ; v ; 

The coefficients of S'^> (m) up to the eight order are given in Table [II]. The intimate coupling 
of this frequency with the sidereal frequency results in the coincidence of the pole in the 
m-series S'J^(m) and S^(m). Numerically Sffl (m Moon ) ~ 0.6035, and the lunar orbit 
sensitivity to the ai perturbation on this frequency is given by |<5 e r| ~ 800 x ai [cm], 
approximately 5.98 times smaller than for the principal sidereal effect. Notice, however, 
that this term in the spectrum of the preferred frame lunar perturbations may bound the 
a± parameter as efficiently as the sidereal term if it turns out that there is significantly less 
noise at this frequency. 



Potential R$ 



A special character of this term is due to its independence on the choice of the gravita- 
tionally preferred frame. It would be theoretically appealing if this term could significantly 
contribute to constraining a\. Unfortunately, we shall demonstrate that the significance 
of this perturbation faces two obstacles: (i) its amplitude is small, and (ii) it acts with a 
synodic frequency, the same as the other phenomena tested through the LLR experiment 
(e.g. the classic equivalence-principle- violation effect; p-|IT||). 

In the context of the Main Lunar Problem we consider a circular solar orbit around the 
Earth. The velocity vo thus becomes — v^ey {vo = n'a') in the rotating (Hill) coordinate 
system introduced in Sec. II. The generating function G reads 



G(w, Dw) = e 2 (-) {[Dw + (1 + m) (1 + w)\ ( - [Dw - (1 + m) (1 + w)} C" 1 } 
with 



£2 — —X21 



2 \ c J ma 
Then the source function can be easily calculated by using ( |3l9] 



(3.24) 



(3.25) 



S (w, w; Dw, Dw) = | Q jc" 1 (1 - m) + ( (1 + m) - 



1 + w 



+ w 



+Dw 



r 1 + c 



1 + w 1 + w 



(3.26) 
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The close similarity with the classical equivalence-principle-violation effect studied in 
consists of the fact that the source function ( |3.26|) excites the odd powers of (, resulting in 
(radial and longitudinal) oscillations of the lunar orbit with the synodic frequency, aliasing 
with the equivalence-principle- violation effect |[l0|Jll|] • 

We have learnt in Refs. [|K],|1]I] that any synodic signal is particularly amplified by the 
presence of a pole singularity occurring for a prograde orbit about 68 % larger than the 
lunar orbit (m cr = 0.19510399 . . .). The corresponding ranging formula reads 



10 



- = -^i 2 1 ) (-)cosr, (3.27) 
a 2m 

where Sj®(m) is a series in the Hill parameter to, whose numerical value for the lunar orbit 



is found to be 1.3022. The amplitude of the synodic oscillation ( 3.27 ) of the lunar orbit thus 



reads \8 e r\ ~ 57 x oc\ [cm], too small to compete with the much greater sensitivity of the 



lunar orbit to the equivalence principle violation term [10,11 



D. Potential R^} 

From Eq. ( |2.4cj ) we see that this perturbing term is equivalent to a variation of the 
gravitational constant G (not to be confused with the generating functions G(w, Dw)) with 
a period of one year (see e.g. [|B|). The analysis of this term involves a small denominator 
A a=m oc — m 2 which, however, cancels out in the radial oscillation |I9"|| . As found in Ref. 0, 
the remaining signal still exhibits an interesting sensitivity to the «i parameter. 

The generating function G equivalent to the reads 

G(w, Dw) = ie 3 (^j (Ce^ - C"V^) , (3.28) 

with 

Iwlfo 

e 3 = ctiC ^ 2 K{m) , (3.29) 



and the resulting source function is 



S K w) = | (") ^- (C m e^ - C m e^) • (3.30) 



Using the previous scheme of solving the perturbation equations we obtain 

5 f r 



(i 



e 3 SS{m)sin[n'(t-t )-<f>] (3.31) 



for the expected contribution to the radial perturbation of the lunar orbit. The first terms 
of the series Sf® (to) are listed in Table |T|. The magnitude of the oscillations (|3.31|) is about 



\6 e r\ ~ 4550 x diC [cm], comparable to the sidereal effect coming from R^- The value 
k (to) Sj® (to) ~ 0.9643 shows that tidal effects are not very important. We learned from 
J.G. Williams (private communication) that the prospects of decorrelating the effect (p. 31 ) 



from other annual effects is high. If this is confirmed, Eqs. ( |3.21| ) and ( |3.31| ) would be the 
two best probes for constraining a±. 



IV. POLARIZATION OF THE LUNAR ORBIT BY A GALACTIC 
DIFFERENTIAL ACCELERATION 

Besides preferred-frame effects, other perturbing forces can be linked to some fixed di- 
rection in space. This would be, in particular, the case if the Earth and the Moon would 
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fall with a different acceleration toward the center of the Galaxy. Years ago this possibility 
has been mentioned in relation with a possible violation of the equivalence principle linked 
to the gravitational binding energy of planets |5B|. More recently, this idea has been re- 



vived within the context of possible strong-gravitational field effects in neutron stars, and 
has led, through the use of existing binary pulsar data, to new tests of strong-field gravity 
| 27j . Still more recently, the idea surfaced again with a different motivation: the possibility 
that the coupling between ordinary (visible) matter and galactic dark matter violate the 
equivalence principle |p8|j . The latter suggestion led to new galactic- related laboratory tests 
of the equivalence principle as well as to a corresponding reanalysis of lunar laser 



ranging data pl|-|33| . Before applying our Hill-Brown algorithm to the latter problem (with 
the result that we find no unexpected amplification), we wish to emphasize that there are 
strong a priori theoretical constraints (using existing observational data) on the conceivable 
magnitude of any "dark matter effect". These constraints diminish, in our opinion, the 
theoretical significance of the results of Refs. ||2l^ , |3T||3"2"|j . 

We assume a field-theoretic framework (as is always the case in recent discussions con- 
cerning possible violations of the equivalence principle; e.g. p8|j30| ). Within such a frame- 
work, any effect on visible matter due to a new (non-Einsteinian) long-range field generated 
by dark matter is necessarily proportional to the product ayai of two coupling constants: 
ay measuring the coupling of the field to visible matter, and otj the coupling to invisi- 
ble (dark) matter. We normalize these coupling constants with respect to the usual Ein- 
steinian coupling so that the effective gravitational constant between bodies A and B reads 
Gab = G*(l ± oiaoib) where is a bare Newtonian constant and where the plus (minus) 
sign holds for a spin (spin 1) mediating field. To fix ideas, let us consider the case of a 
scalar field (our argument goes through in both cases, but is newest in the scalar case, the 
vector case having been already discussed in pi, though with a less stringent constraint on 



ati). Equivalence principle tests probe the composition dependence of the visible couplings: 
a A — 7^ 0. Let us denote by Jab the fractional modification of the average coupling to 
visible matter ay in a differential composition-dependent experiment: a A V ^ — a^P = Jab^v- 
The ordinary tests of the equivalence principle (using visible sources) measure the fractional 
differential acceleration 

/A \ vv 

VZ= — ) = («2° - "STVv = fABOy - (4.1) 



a / ab 

These experiments give us a handle (in practice, an upper limit) on the magnitude of 
ay : ay = (Vab/ fAB) 1 ^ 2 - [For simplicity, we do not put absolute value signs around a, t] 
and /.] Then the tests of the equivalence principle using the same pair of visible bodies, 
and an invisible source (e.g. the galactic dark matter) measure, 

Vab = (—) = (pSP ~ a ( p) a i = fABUyOti . (4.2) 



AB 

Inserting the previous value of ay into ( |4.2| ) yields 

V V A I B = UABY l \ri V A V B) ll ^i. (4.3) 

The main point is, now, that observable facts give not only very stringent limits on i]ab-> but 
also mild limits on aj. Indeed, Damour, Gibbons and Gundlach P4LP5|] have shown, in the 
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case of different scalar couplings to visible and invisible matter, that cosmological data were 
putting the limit = a/2/?/ < 0.71. [See also the later work [36|] which considered only 



the case ay = 0, a>i 7^ 0.] An even more stringent limit comes from gravitational lenses. 
Indeed, in some gravitational lenses one measures three different "gravitational masses": a 
"virial" mass G*(l + aj)M linked to binary interactions, the gravitational mass probed by 
the X-ray-emitting gas G*(l + av«/)M, and the lensing mass G*M (light being uncoupled 
to the new field, be it scalar or vector). The coincidence, within better than 30 %, of these 
three masses in some systems |37j gives the limit ai < (0.30) 1 / 2 = 0.55. [Modulo the sign 
change aj — > —aj, this argument applies to the vector case and therefore improves upon the 
limit c4 vector) < 1 used in pj.] We can apply the above limits to the case of the elemental 



compositions of the Moon (silica) and the Earth (iron core plus silica mantle). A laboratory 
approximation of this case (Si/Al versus Cu) has given 77^ = (5 ± 7) x 10~ 12 |5U| so that, 
at the one sigma level, I^XbI 1 ^ 2 < 3.4 x 10~ 6 . Finally, we get for the maximum Moon-Earth 
(1 — 2, keeping our previous labels) differential acceleration caused by a possible anomalous 
coupling to dark matter (taking into account a further factor m corc /m Eairtil = 0.32) 



VI 



— ) <6.0xlO- 7 v /|/ AB |. (4.4) 
a / 12 

Moreover, the fractional composition-dependence Jab (where A = Si02, B = Fe) is gener- 
ically expected to be small compared to onef]. For instance, in dilaton models |38j one finds 
f A B ^ 1-89 x 10- 5 [(E/M) A -(E/M) B ] where E = Z(Z - l)/(N+ Z) 1 / 3 . Here: Z = atomic 
number, N = neutron number, M = mass in atomic mass units. The only (physically mo- 
tivated) case, we know of, where a a would exhibit a significant fractional variation over the 
periodic table is the case of appreciable coupling to lepton number L, or to B — L = N. 
[In view of the small variation of the baryon to mass ratio (B/M)ab ~ 10~ 3 , any coupling 
involving L with a relative coefficient of order unity leads to essentially the same results.] 
For instance, for a coupling to B — L, one gets /ab = {^N/M) S io 2 ~ {^N/M) Fe ~ —0.076. 
Inserting this figure in Eq. ( |4.4[ ) and using the full galactic acceleration a ~ 1.9 x 10 -8 cm/s 2 , 
one gets 

(Aa)^/ < 3.1 x 10~ 15 cm/s 2 , (4.5) 

which is ten times smaller than the upper limit found in a recent analysis of LLR data 
32 1 . Even if we take \Jab\ ~ 1, we get (Aa)] 7 / ^ 1 x 10~ 14 cm/s 2 , which is three times 



smaller than the result of |[32|| . We conclude that, within what we consider the most natural 
theoretical framework LLR data (and a fortiori laboratory experiments |30|]) do not (yet) 



probe a theoretically very significant domain of values of possible anomalous couplings to 
dark matter. 

Denoting Nc the projection of the unit vector directed toward the galactic center on the 
ecliptic plane and R = (X, Y), the galactic polarization effect is described by the potential 



4 If Iab were larger than one, one should modify our analysis above, and define more carefully the 
average value ay. 
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R G = A G (N G .R) 



(4.6) 



analogous to Eqs. fl2.4p. The parameter A G phenomenologically represents a differential 



acceleration of the Moon and the Earth toward the galactic center. The corresponding 
generating function G is given by 



G = uj 

where 



(1 + w) C 1+m e- i<t>G + (l + w)C (1+m) e"H , (4.7) 



Lb = m 2 4| . (4.8) 
an' 1 

The polar angle 0g = 1.1° gives the angular distance of the galactic center from the lunar 
(and solar) position corresponding to the above chosen new-moon phase at MJD51 521.2. 
Employing Eq. (|3.9[) we obtain the source term of Hill's problem in the following form 



S = _£ e *<fc£-(i+m) _ (4 9) 

Because of the similarity of this function with (pj.l8|) we recover the qualitative conclusions 



of Sec. III.B. The sidereal perturbation of the lunar orbit reads 



a m 

with 



2 ^2 S Um) cos [n {t - t ) - <f> G ] , (4.10) 



r, / x , 75 235 2 127637 , 

o ea i m =1 m H m m 

« 4 384 



4172299 45 
2304 



A more complete set of the coefficients of this series is given in Table 0. The dominant 
m-dependence of this series is again captured by factorizing [^(m)] -1 . The numerical value 
of the series ( |4.11| ) for the lunar orbit is 5 ga i(mMoon) = 0.5050. Clearly, the dark matter 
differential coupling contributes also to the (n — 2 n') -frequency of the radial oscillation of the 
lunar orbit. We do not give here the detailed result, just quoting that its amplitude is about 
5.94 times smaller than the amplitude of the principal galactic polarization contribution 
( |4.1U| ). Finally, the series S ga \(m) shows the same pole, near m cr ~ —0.18407 as the sidereal 
series in ( [3.19] ). 



V. CONCLUSIONS 



The main results of this paper may be summarized as follows: 



We have confirmed, by more detailed computations, previous suggestions |0JT^| that 



LLR data have the potential of constraining the post-Newtonian parameter a% at the 
1 x 10~ 4 level or better. We showed that the preferred frame perturbations associ- 
ated with the «i parameter contribute a large spectrum of frequencies in the radial 
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oscillation of the lunar orbit. The dominant ai-effects occur at frequencies n (sidereal 
effect) and n' (yearly effect) with well-determined phases, and there is a sub-dominant 
effect at frequency n — 2n'. Although the analytical results that we obtained from 
a high-order Hill-Brown algorithm should be accurate enough for fitting purposes, it 
may be advisable to resort to a direct numerical integration of the equations of motion 



see e.g. |19| for the a>i contributions to the equations of motion). 



• We found that retrograde planar orbits with n' / (n — n') = —0.18407 (which have fixed 
perigees in inertial space) exhibit a resonant amplification of preferred-frame effects. 
Putting an artificial satellite near such an orbit could be an efficient (though expensive) 
way of improving the present bounds on ot\. 

• We have extended our analysis to another perturbation linked to a fixed direction in 
space: namely, a possible differential acceleration toward the galactic center. Evi- 
dently, this perturbation exhibits also a pole at m cr = —0.18407 corresponding to an 
orbit "frozen in space". We argue, however, that there are strong a priori theoretical 
constraints on the conceivable magnitude of such an effect. 

A specific suggestion for future work concerns applying our analysis of the singularly 
perturbed spatially frozen orbits to planetary satellites. It is widely known that the solar 
system satellites are submitted to a complicated cosmogonic tidal evolution. It might be 
interesting to study if some of these bodies evolved historically through such a frozen orbit 
configuration yielding indirect limits on the a\ parameter. For instance, one easily verifies 
that several of Jupiter small satellites do lie close to the frozen configuration. A careful 
study of these problems is, however, beyond the scope of this paper. 
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TABLES 



TABLE I. Coefficients s k of the S ( a \' (m) series in powers of m. The percentage pk of the 
contribution of the listed terms to the series for the lunar orbit - m = rriMoon = 0.0808489375 . . . 
- is given in the second column. The last column gives the ratio (sfc_i/sfc). 

k s k p k (sk-i/sk) 






1.000000000000000 






1 


-8.375000000000000 


-67.71 


-0.119402985 


2 


49.37500000000000 


32.27 


-0.169620253 


3 


-268.2473958333333 


-14.18 


-0.184065161 


4 


1444.161892361111 


6.17 


-0.185746070 


5 


-7841.277434172453 


-2.71 


-0.184174314 


6 


42586.57408613037 


1.19 


-0.184125575 


7 


-231375.2814872327 


-0.52 


-0.184058443 


8 


1256984.833107015 


0.23 


-0.184071657 


9 


-6828885.820821555 


-0.10 


-0.184068802 



TABLE II. 


Coefficients s' k of the S'a^ (m) 


series in powers of m. Other parameters 


as in Table I. 


k 


Sk 


Pk 


(Sk-l/Sk) 





1.000000000000000 






1 


-7.166666666666667 


-57.94 


-0.139534884 


2 


40.09305555555556 


26.21 


-0.178750823 


3 


-216.8476851851852 


-11.46 


-0.184890401 


4 


1173.979870756172 


5.02 


-0.184711587 


5 


-6371.734190136318 


-2.20 


-0.184248093 


6 


34622.85594255030 


0.97 


-0.184032600 


7 


-188086.0224857120 


-0.42 


-0.184079899 


8 


1021827.835270500 


0.19 


-0.184068212 
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TABLE III. The same as in Table I but for the (m) series (ratio of the consecutive coeffi- 
cients omitted). 



k 


Sk 


Pk 





1.000000000000000 


— 


1 


-2.000000000000000 


-16.17 


2 


-1.953125000000000 


-1.28 


3 


-5.875000000000000 


-0.31 


4 


12.50295003255209 


0.05 


5 


49.86885579427062 


0.02 


6 


73.34703290903986 


< 


7 


-113.6395324400916 


< 


8 


-670.4626529838962 


< 


9 


-1047.599614057912 


< 





TABLE IV. The same as 


in Table I but for the S gai \(m) series. 




k 


Sk 


Pk 


(Sfc-l/Sfc) 





1.000000000000000 






1 


-9.375000000000000 


-75.80 


-0.106666667 


2 


58.75000000000000 


38.40 


-0.159574468 


3 


-332.3880208333333 


-17.57 


-0.176751256 


4 


1810.893663194444 


7.74 


-0.183549166 


5 


-9846.043167679392 


-3.40 


-0.183920955 


6 


53488.73992091048 


1.49 


-0.184076932 


7 


-290602.0005454238 


-0.66 


-0.184061843 


8 


1578765.334003927 


0.29 


-0.184069155 


9 


-8577020.610670500 


-0.13 


-0.184069201 
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FIGURES 



FIG. 1. Amplitude (in centimeters) C ai of the sidereal oscillation (for ol\ = 1) vs. the 
Hill parameter m (positive for prograde orbits, negative for retrograde orbits). For high orbits 
C ai — ^aeiSai (m)m~ 2 as given in Eq. (3.19). In the case of low orbits we introduce the influence 
of the Earth multipolar structure by adding a quadrupole contribution to the denominator of a 
Pade approximant of Sai(m). Besides the singularly amplified orbit at m cr = —0.18407, two 
celestial bodies are indicated: (i) the Moon (M), and (ii) the LAGEOS satellite (L). 
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